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Abstract 

The near-threshold 12 C(0^~) resonance provides unique possibility for fast helium burning in stars, 
as predicted by Hoyle to explain the observed abundance of elements in the Universe. Properties 
of this resonance are calculated within the framework of the a-cluster model whose two-body and 
three-body effective potentials are tuned to describe the a — a scattering data, the energies of the 
Of and 0^ states, and the 0^-state root-mean-square radius. The extremely small width of the 
0^ state, the 0^~ — > Of monopole transition matrix element, and transition radius are found in 
remarkable agreement with the experimental data. The 0^-state structure is described as a system 
of three a-particles oscillating between the ground-state-like configuration and the elongated chain 
configuration whose probability exceeds 0.9. 

PACS numbers: 21.45.-v, 21.60.Gx, 21.10.Tg, 24.30.Gd, 27.20.+n 
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The extremely narrow 0(0^) resonance just above the 3a threshold was predicted by 
Hoyle l| to explain the observed abundance of elements in the Universe. Following the 
prediction, the near threshold 3a ( 12 C(02~)) and 2a ( 8 Be) resonances are crucial for suffi- 
ciently fast carbon production in stellar nucleosynthesis which goes through the reactions 



3a -> 8 Be + a ->■ 12 C(0 



— >• 
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C + 7. The experimental confirmation came soon |2j and 



properties of this state have been fairly well studied experimentally during the following 
50 years. In particular, the extremely small width F and the — > Of transition density 
(the monopole transition matrix element M 12 and the transition radius R12) were deter- 
mined. Besides the resonance triple-a reaction, a consistent treatment of the three-body 
near-threshold dynamics is ultimately necessary to study the non-resonance reaction, which 
is of importance at low temperature and high helium density, as occurs in accretion on white 
dwarfs and neutron stars [3j. As the experimental data at these conditions are lacking, the 
reliable theoretical results are especially desired; nevertheless, the calculated non-resonance 
reaction rates Q-Q differ by many orders of magnitude. There is vast literature on the 12 C 
states above the 3a threshold, and it is possible to notice just some of the recent papers, 
e. g., with concern of astrophysical asp ects [z-9], distribution of decaying particles 10, U|, 
dilute-gas variational calculations [l2j, and different approaches to precise description of 



the 12 C properties 13Nl7l|. Besides much interest in the 12 C(0^) state owing to its role in 



the astrophysical applications, the structure of this state has been widely discussed for a 
long time, in particular, a linear chain configuration was suggested in 18] and a triangle 
configuration was considered in [lj3]. 

The following reasons impede the theoretical description of the 12 C(0^~) resonance and, 
more generally, of the 3a low-energy processes. Application of the twelve-nucleon calcula- 
tions is hardly feasible for the unstable near-threshold 0^ state; the extremely difficult task 
is to describe the a-cluster structure and complicated three-body dynamics of a-particles. 
On the other hand, although the a-cluster model is a natural alternative and a number of 
works have elaborated in this direction, it remains to prove that the model would provide 
the desired accuracy. One of the crucial difficulties stems from the three-body decay of the 
0^ resonance, which requires the construction of the asymptotic wave function for three 
charged particles. 

The present Letter is aimed at reliable calculation of the 12 C(02~) resonance within the 
framework of the a-cluster model with particular intention to obtain its extremely small 
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width T and the experimentally observable monopole transition matrix element M\i and 
transition radius Rn- Moreover, the reliable a-cluster calculation will be used to elucidate 
the O2 -state structure. The a-cluster model is applied also to calculate the 0+ state de- 
spite that its cluster structure is not justified. First, this calculation provides a consistent 
treatment of the 0^-0^" transition as the wave functions of both states should be calculated 
within the same approach. Next, it is of importance for a correct description of those effects 
which originate from the a-particles structure. An additional goal of the 12 C(02~) calculation 
is to construct the effective potentials of the a-cluster model which is known to be successful 
in different nuclear problems. The focus of the present calculation is the precise solution of 
the three-body problem which allows one to show an ultimate ability of the a-cluster model 
in description of the fine 12 C characteristics. 

For a precise description of the 12 C nucleus, firstly, the effective three-body interaction 
is introduced to take into account the finite size and nucleonic structure of the a-particles 
including the effects of wave- function antisymmetrization under nucleons. Next, the as- 
sumption on the "frozen" nucleon distribution in the a-particles is used to calculate the 
experimentally observable nucleonic structural characteristics, e. g., the root-mean-square 
(r. m. s.) radii and the transition radius Ryi- This extended a-cluster model leads generally 
to a satisfactory description of the 12 C nucleus 
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T and M12 
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20| . In particular, the calculated 



14[ are in reasonable agreement with experimental data (up to a factor two) 



even for the simplified local a - a interactions. 

The effective potentials of the a-cluster model are constructed under the following nat- 
ural requirements. The two-body potential is chosen to reproduce the low-energy elastic- 
scattering phase shifts in the lowest even-parity waves and to fix at the experimental values 
both energy and width of the narrow a-a resonance (the 8 Be ground state). The latter con- 
dition is essential for a correct description of the multi-dimensional potential barrier whose 
penetrability determines the extremely small width T of the 12 C(0 2 f ) state. Furthermore, 
for any reasonable calculation of the width T, the O^-state energy should be precisely fixed 
to the experimental value by a proper choice of the effective three-body potential. As the 
ground-state wave function is needed to treat the 0^-0^ transition, the effective three-body 
potentials should provide also the experimental values of the energy and r. m. s. radius of 
the Of state. 

The Schrodinger equation reduces to the system of coupled hyperradial equations (HRE) 
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for the channel functions f n {p) by using the expansion of the three-body wave function in a 
set of eigenfunctions $„(«, 9; p) on a hypersphere at fixed p, 



oc 



y = p-* l2 Y,Up)®n{^d ]P ). (1) 



n=l 



Here the hyperspherical variables 0<p<oo,0<aj<7r, and < 9i < 7r are defined by 
= pcos(aj/2), ?/j = psin(o;j/2), and cos#j = (xj ■ yi)/(a?ij/i), the scaled Jacobi coordinates 
are x, = — r^, = (2rj — r,,- — rfc)/\/3 and is the position of the ith a-particle. Both 
the total wave function \l/ and the eigenfunctions & n (a,6;p) are symmetrical under any 
permutation of the a-particles. In more detail, the method of calculation and the numerical 
procedure are presented in 3, Q . 

The two-body interaction includes the Coulomb part, 4e 2 /x, and the short-range part of 
the simple Ali-Bodmer form, 

V(x) = V r e~ x2/ ^ + V a e~ x2/ ^ , (2) 

in the s, d, and g partial waves. Following the above discussed requirements, the parameters 
of the two-body potentials are chosen to fit the low-energy a - a elastic-scattering phase shifts 
(up to E cm = 12 MeV) and to reproduce the experimental energy E2 a = 92.04 ± 0.05 keV 



and width 7 = 5.57 ± 0.25 eV of the 8 Be ground state [2l|, |22J. These conditions do not 
uniquely fix the two-body interactions and to elucidate its role two variants are considered 
for each of the s-, d-, and g-wave potentials whose parameters are listed in Table [H Given 
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TABLE I: Parameters V r ^ a (MeV) and [i r ^ a (fm) of the a-a potential ([2]) for two variants in each 
of the / = 0, 2, 4 partial waves. 

the two-body potential, the corresponding three-body potential of the Woods-Saxon form, 

V 3 (p) = V [l + e^ b ]~\ (3) 
4 



is chosen to fix the ^ -state energy, the O2 -state energy, and the 0+ -state r. m. s. radius at the 



23|, and i?i = 2.48 ± 0.002 



experimental values E x = -7.2747 MeV, E 2 = 0.3795 MeV 
251 ]. Five of eight possible combinations of the two-body potentials under consideration, 
which provide the desired Ei, E 2 , and Ri, and fitted parameters Vq, a, and b of the three- 
body potential (jSJ) are presented in Table [TT1 
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TABLE II: The 0^ state width T (eV), the r. m. s. radius R 2 (fm), the monopole transition matrix 
element M\ 2 (fm 2 ), and the transition radius R\ 2 (fm). Parameters of the two-body potential 
indicated according to Table U and parameters of the corresponding three-body potential ([3]) are 
V (MeV), a (fm), and b (fm). 

A set of eigenfunctions Q n (a,6;p) and all the terms in the HRE are obtained by using 
the variational method with a flexible basis of symmetrical (under any permutation of three 
particles) trial functions. For small p < 20 fm, sufficient precision is achieved on the basis 
containing 192 (corresponding to K max = 90) symmetrical hyperspherical harmonics (SHH). 
With increasing p, a number of SHH providing the desired accuracy enormously grows due 
to the a + 8 Be cluster structure of the wave function. Therefore, the desired precision is 
obtained with the basis containing a set of 108 SHH (corresponding to K max = 66) and four 
additional trial functions (fj = J2i e W {~Pj x D- 

Whereas the simple zero asymptotic boundary conditions in the HRE are used in the 
calculation of the 0^ state, the complicated scattering problem for three charged particles 
should be solved to describe the 0^ resonance. Fortunately, the a + 8 Be two-cluster struc- 
ture of the wave function at large, though finite, distances is corroborated in the present 
calculations, which allows one to circumvent the tremendous calculation of the three-charged- 
particle wave function. The a + 8 Be structure follows from the form of $i(a, 8; p) which is 
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close to the symmetrized combination J2i 4>{xi) of the 8 Be wave function 4>{x) in a wide region 
20 fm < p < 70 fm, where the overlap integral on a hypersphere Yli 0(^0) exceeds 0.977. 
This form of 8; p) is demonstrated in Fig. [TJfor p = 30 fm, where three peaks represent 
the functions <f)(xi). The two-cluster structure is approved also by the form of the first- 




FIG. 1: The eigenfunction $i(a, 6; p) at p = 30 fm. 



channel effective potential in the HRE, which to good accuracy is a sum of the 8 Be energy E 2a 
and the Coulomb interaction between a and 8 Be, 2t]e 2 /p in the interval 40 fm < p < 70 fm. 
Indeed, fit to this dependence gives 89.29 keV < E 2a < 90.4 keV and 6.73 < r] < 6.76 which 
are close to the expected values E 2a = 92.04 ± 0.05 keV and i] = 8/y/3 ~ 6.65. Thus, the 
0^"-state properties are calculated by using the two-cluster (a + 8 Be) asymptotic boundary 
condition for the first-channel function f\(p) at a large, though finite, hyperradius. More 
precisely, the asymptotic form is fi{p) oc F (r]/kr b , kp) + tan 5(E)G (ri / kr b , kp), where 5(E) 
is the scattering phase shift, F and G are the Coulomb functions, k 2 = (m a /h 2 )(E — E 2a ), 
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and r b = h 2 /(m a e 2 ). The O^-state energy E 2 and width T are obtained by fitting 5(E) to 
the Wigner energy dependence. 

The assumption on the fixed nucleon distribution in the a-particle leads to the representa- 
tion of the twelve-nucleon wave function as a product of the three-body wave function \l/ and 
the internal wave functions of the three a-particles. Using this form, one obtains from the 



usual definitions 
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26| that the monopole transition matrix element is expressed merely 



via the three-body wave functions, M\ 2 = \p 2 \ ^^), whereas the r. m. s. radii and the 
transition radius depend also on the a-particle r. m. s. radius R a = 1.681 ± 0.004 fm I2T 



viz., R 2 = \p 2 \ *W) + Rl and R 
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( " (1)|p4(3 ~;S: m2e)l " (2)> + f K ^ calculation 
of M12, R 2 , and R\ 2 the ultra-narrow O^-resonance is treated as a bound state, whose wave 
function is taken as a scattering solution at the resonance energy E 2 normalized on the 
interval p < pt, where pt ~ 45 fm is the typical extention of the potential barrier. 

The solution of the five-channel system of HREs provides the desired accuracy for T, M 12 , 



R12 , and R 2 which are presented for all the selected potentials in Table [TT1 Amazing 
T and M\2 coincide with the measurements 8.5 ± 1.0 eV 



23| and 4.396 ± 0.27 fm 2 [26j within 



y, both 



ightly overestimates the experimental 
of significantly improved accuracy 



the experimental errors, while the calculated R\ 2 s 
value 4.396 ± 0.27 fm 26|. The recent experiment 
gives M12 = 5.47±0.09 fm 2 and R\ 2 = 4.59±0.16 fm, still the calculated M\ 2 is in agreement 
with these highly accurate data, whereas R\ 2 agrees even better. Thus, the quite natural 
requirements on the simple two-body and three-body effective potentials of the a-cluster 
model turn out to be sufficient for a precise description of the lowest 12 C(0 + ) states. The 
reliable value of the O^-state r. m. s. radius is found with a small variance about R 2 = 3.57 fm 
which exceeds R\ more than 1 fm. Note that the exact experimental data on R 2 are not 



no 



er diffraction in the 



available and the efforts to extract its value, e. g., from the Fraun 
inelastic a + 12 C scattering, are affected by the model assumptions [29 ]. 

A stark success in the calculation of the T, M\ 2 , and R\ 2 implies that the a-cluster 
model provides a good description of the O2 state's structure. The structural properties 
become more conspicuous, if represented in the variables £ and fi defined by the expressions 
sin£ = sinaisin^j and cos£ cosipi = cosaj. As the inertia momenta l\ = \m a p 2 cos 2 (£/2), 
I 2 = \m a p 2 sin 2 (£/2), and J 3 = \m a p 2 are related to p and £, the probability distribution of 
ct-particles P(p, £) = ^p 5 sin2£J dip\^/(p, £, <^)| 2 is plotted in Fig. [2]to represent the features 
of the O2" state. The sufficiently wide main peak corresponds to different configurations of 
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FIG. 2: Probability distribution of a-particles P(p,£) for the 12 C(02 ) state. 

large size, mostly satisfying the condition J 2 3> J 1; in particular, P(p, £) takes the maximum 
for the ratio ~ 18 and p = 6.3 fm. In other words, the elongated chain configuration 

is the most probable shape of the state. The peak extends to a ridge in the direction 
of increasing I2/I1 which clearly represents the a + 8 Be decay mode. The minor peak 
is quite similar to the 0^-state probability distribution whose form is consistent with the 
equilateral triangle configuration. The two peaks are separated by a ravine pointing out 
at the position of the nodal surface which separates the regions of different sign of the 
wave-function. This means that three a-particles in the state oscillate between the 
ground-state-like configuration and the linear chain configuration which prevails in the total 
probability, whereas the weight of the minor-peak region is about 7-9%. A detailed and 
lengthy discusssion of the state structure is beyond the scope of this Letter and will be 
given elsewhere. 



S 



In summary, it is shown that the a-cluster model, in spite of its relative simplicity, 
provides a consistent description of the lowest 12 C(0 + ) states. Even for the simple effective 
potentials describing the a-at scattering data, the energy and width of 8 Be, the energies of 
12 C(0^ 2 ) states, and the 0^ r. m. s. radius, the calculated state width T, the monopole 
transition matrix element M 12 , and the transition radius R\2 are in excellent agreement 
with the experimental data. The wave function of the state corresponds to the system 
of three a-particles which oscillate between the equilateral triangle and the linear chain 
configurations, the latter contributing to the total probability about 90%. The three-body 
calculation directly approves the sequential decay mechanism ( 12 C — > a + 8 Be — > 3a) of the 
0^ state in agreement with the experiment 30]. 

The present results indicate that the theoretical calculations within the framework of the 
a-cluster model are accurate enough to keep pace with the experimental data. There is 
enough room for further refinement of the a-cluster model which could be used to describe 
both the improved experimental data and those reactions which depend sensitively on the 
fine details of the wave function. Important examples could be the a-a bremsstrahlung 
and the (a, a) reactions. Furthermore, the present approach is promising to study the 
triple-a reaction (3a — > 12 C) at low energy below the three-body resonance that provides 
an opportunity for unified treatment of the crossover from the resonant to the non-resonant 
mechanism of the reaction. 
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